ELEMENTARY MATRIX DECOMPOSITION AND THE COMPUTATION OF 
DARMON POINTS WITH HIGHER CONDUCTOR 



XAVIER GUITART AND MARC MASDEU 

Abstract. We extend the algorithm of IDG02I and IDP06I for computing p-adic Darmon points 
on elliptic curves to the case of composite conductor. We also extend the algorithm of IDL03I for 
computing ATR Darmon points to treat curves of nontrivial conductor. Both cases involve an 
algorithmic decomposition into elementary matrices in congruence subgroups Ti(?t) for ideals 91 
in certain rings of 5-integers. We use these extensions to provide additional evidence in support 
of the conjectures on the rationality of Darmon points. 



1. Introduction 

Let E be an elliptic curve over Q of conductor pM, with p a prime not dividing M. Let K be 
a real quadratic field in which p is inert and all the primes dividing M are split, and denote by 
Hp = ¥ 1 (K P ) \ P 1 (Q P ) the i\~ p -points of the p-adic upper half plane. 

A construction of Darmon |Dar01| associates to every t <E T~L p a local point P T E E(K p ), which 
is defined as a certain period of the modular form / corresponding to E under the Modularity 
Theorem. The points P T are conjectured to be rational over ring class fields of K , and to behave 
in many aspects as Heegner points. 

An algorithm for the effective calculation of p-adic Darmon points was given in |DG02j and 
was improved in jDP06j . Some P r 's were computed in concrete examples and checked to be p- 
adically close to global points, providing extensive numerical evidence in support of the conjectures. 
However, due to the restrictions imposed by the algorithm, only elliptic curves of prime conductor 
p (that is, with M = 1) could be treated. 

In the articles }DG02] and |DP06j it is crucial to assume that M = 1 when applying the "contin- 
ued fraction trick" (see |DG02( p. 42]) in order to transform certain semi- indefinite integrals into 
double integrals. The present article provides a different procedure for performing this step when 
M > 1 provided that E satisfies a mild condition (see Assumption 13.21 for its precise statement). 

In other words, we extend the algorithm of |DG02| and jDP06j to a much larger class of curves. 
As an application, we compute p-adic Darmon points on curves of composite conductor and we 
check that they are close to global points, which provides new experimental evidence in support of 
the validity of Darmon's construction in its stated generality. 

The method for transforming semi-indefinite integrals into double integrals is based on an al- 
gorithmic decomposition into elementary matrices in congruence subgroups Ti (91) for ideals 9T in 
certain rings of S-integers. This can be seen as an effective version of the congruence subgroup 
problem. In particular wc improve on a theorem of Cooke- Weinberger (see Corollary 12. 4|) . 

In addition, essentially the same method of elementary matrix decompositions can also be applied 
to ATR points, a different instance of Darmon points introduced in |Dar04i §8] and |DL031 §4] for 
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elliptic curves over totally real fields. Although the construction of ATR points is radically different 
(for instance, they are defined by means of complex periods), their explicit computation has some 
formal similarities with respect to the p-adic setting. In particular the methods used until the 
present have all used a "continued fraction trick" which in this case only applies to curves of trivial 
conductor (cf. jGarlll p. 108] for a discussion of this issue). The present article provides also a 
method for computing ATR Darmon points in curves of non-trivial conductor. 

The rest of the article is organized as follows. In Section [2] we introduce an algorithm for 
computing elementary matrix decompositions in certain congruence subgroups. We state it in a 
level of generality so that it can be applied both to the p-adic and the ATR setting. In Section [3] 
we recall the definition of p-adic Darmon points, we discuss the algorithm for computing them in 
curves of composite level and we include some tables of numerical computations performed using it. 
Finally, in Section [4] we briefly recall ATR points and we make explicit the method for computing 
them in curves of non-trivial conductor, as well as a detailed example of a numerical verification of 
Darmon's conjecture in this case. 

Acknowledgments. It is a pleasure to thank John Voight for helpful comments on elementary 
matrix decompositions, Robert Pollack for clarifying to us some details on his implementation of 
the overconvcrgent modular symbols algorithm, and Henri Darmon for many valuable observations. 
Guitart wants to thank the Max Planck Institute for Mathematics for their hospitality and financial 
support during his stay at the Institute, where part of the present work has been carried out. 



2. Preliminaries: elementary matrix decomposition in Pi 

Let F be a number field with ring of integers O, and let S be a finite set of places of F containing 
the archimedean ones. Let Os denote the subring of F consisting of those elements whose valuation 
is non-negative at all the places outside S. 

For an ideal 9t of Os, we denote by Pi(0 f l) the subgroup of SL 2 (0s) defined as 

(2.1) r 1 (9T) = | 7 eSL 2 (O s ): 7 = (I (mod 9?)}, 

and by the subgroup of Ti(yV) generated by the elementary matrices of the form 

(2.2) Q fj with y <G 9t and ^ with x G Os- 

If the group of units Og is infinite, then Pi(9t) = Ei^; see for instance |Vas72j for a proof of 
this result, and also for its relation with the Congruence Subgroup Problem. Therefore, when Og 
is infinite every matrix in ri(9t) can be factored into a product of elementary matrices of type 
(|2.2|) . However, the proof given in |Vas72j is not explicit, so it does not give rise to an algorithm 
for systematically performing such decomposition. In this section we see that the results and 
techniques introduced in |CW75j can be adapted to provide -assuming GRH- such an algorithm, 
in the particular case where F has at least a real archimedean place. Thus from now on we will 
assume that Og is infinite and that F has at least a real place. 

The following lemma is stated and proved in |BMS67| Lemma 2.2 (b)] using the notation of 
Mennicke symbols. For the convenience of the reader we restate it in terms of explicit matrix 
formulas. 
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Lemma 2.1. Le^7=("^) be an element in Ti^). Suppose that c = u + ta for some unit u £ Og 
and some t £ Os, and let T be the matrix 

( 1 o\A u-i\( 1 OUl 
\u(a-l) l) \0 1 J\t(a-1) l) \-c 1 

Then T7 = (01) f or some x £ Os- In particular, 

c + t(l-a) lMo 1 JU(l-a) 1J \0 1 



(2.3) 7 = 



Proof. A direct computation shows that the first column of T7 is ( J ) . Since T7 belongs to ri(9T), 
we see that T-y = (Jf) for some x £ Os- □ 



Observe that, since a — 1 and c belong to 91 , identity (|2.3[) already expresses 7 as a product of 
elementary matrices of type (|2.2[) . The next step is to show that, assuming GRH, one can reduce 
to the case where a is congruent to a unit mod c by multiplying 7 by an elementary matrix. Before 
stating the result of |CW75j that grants this, we recall some terminology related to ray class groups. 

Let m be an ideal of O, with factorization into prime ideals of the form 



m=nv 



p|m 

Let J m be the multiplicative group of ideals of O relatively prime to m, and let 

i*i,m = {x £ F : v p (x - 1) > m(p) for all p | m}. 
There is a natural map i: -Fi, m — > I m given by x > (x). The quotient 

C m = I m /i(F hm ) 

is the ray class group modulo m. The following result is |CW75l Theorem 2.13]. 

Theorem 2.2 (Cooke- Weinberger). Assume GRH. Let m be an ideal in O, and let a be an ideal 
class in C m . Then the set of prime ideals q contained in a such that the reduction map 



0* — ► (Os/qOs 

is surjective has positive density. 

From this we obtain the main result of this section. 

Theorem 2.3. Let 7 = (" d) oe an element o/r^lt). Assuming GRH, the following algorithm 
terminates and computes an expression of 7 as a product of elementary matrices. 

(1) Iterate over the elements A £ O to find A such that a' = a + Ac generates a prime ideal and 
(2.4) Os — > (Os/a'Os) X is surjective. 

(2) 5etV = (SJ)7, and let Y = («,'*',). 

(3) Iterate over the elements u £ 0$ until finding u such that 

c = u (mod a). 

(4) Use Lemma \2~7\ to find an expression of 7' as a product of elementary matrices 
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Proof. We can, and do, choose e £ Og such that ea and ec belong to O. By applying Theorem l2.2l 
with m = (ec), we see that there exists a prime ideal q in the same class of [(ea)] £ C m such that 
Og — » (0/a'0s) x is surjective. Since (ea) is integral and principal, we see that q is also principal. 
Therefore, q = (g) for some q = ea (mod (ec)), so that q = ea + Xec for some A £ O. But q and 
a' = a + Ac generate the same ideal in Os, so that Og —> (0/a'Os) X is surjective. This justifies 
that step |JlJ) of the algorithm can be accomplished. 

Now the class of c' in (Os/a') x can be represented by some unit u £ Og , and this justifies step 
©. □ 

In the particular case 7t = Os, the result in |CW75[ Th. 2.14] asserts that any matrix 7 £ 
SL2(0s) can be expressed as a product of at most 7 elementary matrices. Theorem 1 2 . 31 above and 
expression (|2.3j) give the following generalization to arbitrary ideals 91, which also slightly improves 
the number of elementary matrices needed to 5. 

Corollary 2.4. Under the assumption of GRH, every matrix in T\(^!t) is a product of at most 5 
elementary matrices of type (|2.2|) . 



3. Computation of p-adic Darmon points on curves with composite conductor 

In this section we explain the algorithm for computing p-adic Darmon points in curves of compos- 
ite conductor and some of the related computational issues. After briefly reviewing the definition 
of p-adic Darmon points in £13.11 in £13.21 we describe the algorithm to transform the semi-indefinite 
integrals appearing in the definition of these points into double multiplicative integrals. In £13.31 we 
explain an efficient way of computing these double integrals, and finally in £|3.4l we comment on the 
calculations that have been carried out in support of Darmon's conjecture. 

3.1. Review of p-adic Darmon points. Our presentation of the necessary background and the 
definition of p-adic Stark-Heegner points follows closely [DP061 §1], to which we refer the reader 
for more details and, in fact, for an excellent account of this material in the prime level case. 

Let E be an elliptic curve over Q of conductor N = pM, with p a prime not dividing M. Let 
f( z ) = En>i flne 27rmz be the weight two newform on Fq(N) whose L-series coincides with that of 
E. The coefficient a p is 1 (rcsp. —1) if E has split (resp. non-split) multiplicative reduction at p. 

Let R be the order in M2(Z[l/p]) consisting of matrices that are upper triangular modulo M, 
and let T = denote its group of units of determinant 1. Let if be a real quadratic field in which 
p is inert and all primes dividing M are split. Set Hp = P 1 (K p )\¥ 1 (Q p ), in which T acts by Mobius 
transformations. The p-adic Darmon point construction yields a map 

r\n p — > e(k p ) 

T .— > P T! 

given in terms of certain p-adic integrals, and whose definition ultimately relies on the Z-modular 
symbol attached to E. 

3.1.1. Measures attached to modular symbols. If V is a Z-module, a V -valued modular symbol tp is 
a map 

ip : P 1 (Q) x P 1 (Q) — ► V, (r, s) h-> ip{r -> s} 

such that 

ip{r -> s} + (p{s -^t} = ip{r -> i) for all r,s,t£ ¥ 1 (Q). 
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For e {±1} we denote by If : P X (Q) x F 1 {Q) -> Z the Z- valued modular symbol attached to 
E and u^. That is to say 



I f {r -> s} = 



J r s Rew/ if Woo = +1) 
f^lmuf if Woo = -1, 



where = 2irif(z)dz and f2 + , f2 _ g M >0 are the unique periods with the property that the map 
// defined in this way takes values in Z and in no proper ideal of Z. To simplify the exposition 
we assume for the rest of the article that Woo = 1, but the construction works very similarly for 

Woo = -I- 

A 7L-valued measure on P 1 (Q P ) is a finitely additive function /i from the set of compact open 
subsets of P 1 (Q P ) to Z. In [DarOlj Darmon attaches to each pair r,s £ P 1 (Q) a Z- valued measure 
H f {r -> s} on P (Qp) with total measure by defining 

H f {r -> s}(-fZ p ) := Ifi^r -> 7 ^ 1 s}. 

This is enough to define fJ.f{r — »• s} for all compact open U C Pi(Qp), because either U or F 1 (Q P )\U 
are of the form 7 Z P for some 7 £ T. 



3.1.2. Double multiplicative integrals. If /i is a continuous function on Pi(Q p ) and fj, is a measure 
on P 1 (Q P ) then the integral J* P i/q j h(x)d/j,(x) is defined by means of the Riemann sum 

/ h(x)dfj,(x) = lim > h(x a )/i(U a ), 

M{Q P ) U={U a }^ 

where the limit is taken over increasingly finer coverings U of P 1 (Q p ) by compact open subsets U a , 
and x a is any point in U a . If \i is Z- valued the multiplicative integral can be defined by replacing 
the Riemann sum by a Riemann product: 

/ h(x)du(x)= lim TT^(a:a) MC/ot) - 

JP(Q P ) a 

For ri,T2 € Hp and r,s£ P 1 (Q P ) Darmon defines a JsT p - valued double integral as 

/ f UJ f = f log ( ~ — ~ ) d M r ^ s ^ x ^ 

where log denotes a fixed branch of the p-adic logarithm. Since Hf{r — > s} takes values in Z a 
K* -valued double multiplicative integral can be defined as 

/ [ UJ f = I ( - — — ] dn f {r -> s}(x). 

Jn Jr JP 1 (Q P ) V X_T 1/ 

These two integrals are related by the formula 
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3.1.3. Semi-indefinite integrals. The double multiplicative integral satisfies the usual additivity 
properties with respect to the limits, as well as the T-invariance property 

(•JT2 /-7 s r T ' 2 



r7T 2 l-ys r-T-2 rS 

x / ujf = i / ujf for all 7 e T. 

J 7Ti J 7r Jti J r 



'7T1 J 7r J ri 

Therefore, it gives rise to a group homomorphism 

Int* : (Div°(H p ) ® Div°(P 1 (Q))) r — ► X p x , 

where the subscript T denotes the subgroup of T-coinvariants. Let Z[r] denote the group ring of T 
and let Zr denote the augmentation ideal, defined by the exact sequence 

— > I T — > Z[T] — > Z — ► 0. 

Tcnsoring with Zp over Z and taking T-coinvariants gives 

— > ZT r — > (J r ® J r )r ^ (Z[T] ® 7 r )r — > (Jr)r — > 0. 

Since Zp is generated over Z by elements of the form 7 — 1, choosing base points r € Hp and 
iff 1 (Q) one can define integration maps 

Int* x : (Z r ® Zr)r —j. Zf* 

determined by 

Int x x (( 7o - 1) ® ( 7l - 1)) = / / w/, for 70, 7l G I\ 

Letting A = I TjX (Kp) C ifp yields a well denned map 

(3.1) Int*,: Im(r)— ^/A. 

The group (Zp)p ~ r a (, is finite, say of exponent ep. If y = jx € P 1 (Q) is in the same T-orbit as x, 
define 

rv 



e r ujf := Int* x (e r ■ 1 ® (7- 1)) G A'* /A. 
This semi-indefinite integral satisfies the following properties: 

/T PS pT ft PT pt 

/ epw/ x i / erujf = f / epw/, for all r € ?Z P , r,s,t £ Tx; 

/ er^f * / erw/ = x / erw/, for all n, r 2 £ r, s G Tx; 

/*7 r /*7 S /* r r s 

(3) x / erWf = f / epw/, for all 7 G T. 



3.1.4. Darmon points. Let q denote the p-adic period of E and let $Tato: K^/q — > E(K p ) be 
Tate's uniformization map. 

Given r G % p , let C T be the ring of matrices in R that have the vector (r, 1) as eigenvector, 
which is isomorphic to a Z[l/p] order of A. Let H T denote its ring class field which we can, and 
do, view as a subficld of K p by choosing a prime of ZZ T above p. The stabilizer T T of r in T is a 
cyclic group of infinite order isomorphic to O* 1 / (±1). Let "f T be a generator of IV, and define 

/T P'YtOO 

Conjecture 3.1 (Darmon). The local point P T = $Tatc(Z T ) belongs to E(H T ). 
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3.2. Computation of semi-indefinite integrals. In order to effectively compute the points J T 
one needs to compute the semi-indefinite integrals 



The method used in |DG02j and |DP06| boils down to using properties (Q]), © and ([3]) of semi- 
indefinite integrals to express them in terms of double integrals, which can be effectively computed 
either via Riemann products as in |DG02j or, more efficiently, via overconvergent modular symbols 
as in |DP06j . 

The algorithm for expressing semi-indefinite integrals in terms of double integrals of }DG02] and 
|DP06| is based on the continued fraction algorithm, and it only works under the assumption that 
AI = 1 (i.e., that E has conductor equal to p). In this section we introduce an algorithm that, 
assuming GRH and under some mild assumptions on E, works for all levels M. In Section [3~3l we 
will sec that the resulting definite double integrals obtained by this method can also be computed 
using overconvergent modular symbols, by suitably adapting the techniques of |DP06j . 

From now on we assume that the modular form / satisfies the following condition, which is 
non-vacuous only when M is prime. 

Assumption 3.2. There is a d > 1, d \ M such that / has eigenvalue 1 with respect to the 
Atkin-Lchner operator Wd- 

Remark 3.3. Under Assumption 13.21 the double multiplicative integral is also invariant under the 
matrix w d = ( _° d J ) : 



Let f be the subgroup of PGL2(Q) generated by T and Wd- Then, by replacing T by f in the 
argument of Section 13.31 one can extend the definition of the semi-indefinite integrals § J er^ / to 
all pairs r,s e P 1 (Q) lying in the same T-orbit. 

Let Li be the congruence subgroup of SL2(Z[l/p]) defined as 



Therefore, replacing P T by a multiple of it if necessary we can always assume that j T belongs to Ti 
(but see also Remark l3"^ij) . 

Observe that Ti is one of the groups treated in Section [2] Indeed, with the notation as in that 
section, if we let F = Q, S = {oo,p} and *Tt = M ■ Z[l/p] we have that Ti = ri(9T). In particular, 
the algorithm described in Theorem 12.31 gives an effective method (under our running assumption 
of GRH) for computing a decomposition of 7 r of the form 



(3.2) 





F 1 ={ 7 er: 7 =(it) (modM)}cr. 
Using the properties of the multiplicative integral it is easy to see that 
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In particular, Li and Ui belong to V . Then, for G € T we have that 

/r pUiGoo pU i 1 r pGoo 

I e T u} f = j I e r uj f 
J QC J J DO 

and 

LiGoo pT pO ft pLiG-oo 



(3.5) 



loo J Jo 

r L i 1 ' r fGoo 

erUJf x x / us+ 

oo J JO 



fi 



rLT -t roo fL~ -r r Goo 



erujf x x / erujf x x ; er^>f 

oo J JO J J oo 

oo pL7~ -r rG-oo 

erUJf xx / erujf- 

o J Joo 

In view of decomposition (|3.3p , repeated application of (|3.4j) and (|3.5[) transforms the semi- indefinite 
integral (|3.2[) into a product of double multiplicative integrals. Observe that Assumption 13.21 is 
crucial in (|3.5j) because of the integral f T f erw/. Indeed, the cusps and oo arc never in the 
same T-orbit if M > 1, but u^oo = so they are in the same T-orbit and thanks to Remark 13.31 
the integral f T erw/ is well defined. 

Remark 3.4. As we already mentioned, one can overcome the fact that 7 T = ( ™ b d ) does not generally 

belong to Ti by computing an appropriate power {J T ) m . In some cases one can apply an alternative 

procedure instead, which turns out to be more convenient in the actual computations and it allows 

for the computation of J T itself. Namely, if a = p n (mod M) for some integer n, then the matrix 

(v~ n 0\ 
g = n n I belongs to T and 



p n 



with gj T belonging to Ti 



ffTOO pgr rgj-roo 

J T = i / erevf = X / er^Jf 



3.3. Computation of the definite double integrals. As we have seen in Subsection 13.21 the 
computation of the period J T is reduced to products of integrals of the form 

(3.6) / (^Adu f {r^s}(x). 

Write n for the measure Hf{r — > s}, and consider a decomposition of P X (Q P ) into a disjoint union 
of L open balls of the form 

(3.7) F 1 (® p ) = \Jg i -Z p , 9i =h b j)eGL 2 (Q), 
which will be fixed later. This yields in tnrn a decomposition 
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Fix such an i and let g = gi be the corresponding matrix, written g = ( " ^ ) ■ We are thus reduced 
to calculating 

' X — T2 



n 



dfj,(x) 



Apply the change of variables x = g ■ t to get 

\x-nj Jz p \g-t-nJ 

Let log p be the unique homomorphism K p — > K p such that log p (l — t) = —J2^Li^ n / n an d 
log p (p) = 0. It is surjective, with kernel 

ker (logp : K* -> K p ) = p z X U, 

where U is the group of roots of unity in K p . Suppose that we can express the integrand as a 
power series in t of the form 

0.8) (^)=«o(i + E;^)> 

with OL n belonging to Z p for all n > 1. Then the expression in p.8j) converges for i g Z p and is 
constant modulo pM^oH 1 . Therefore the expression of (|3.6p can be determined modulo this power 
of p by performing L evaluations. The logarithm log ( J T ) of the period is evaluated by noting that 

\og p (a (l + £ a„p"n) = log P ao-E ~ ~ • 

n— 1 

Interchanging the infinite sum with the integral, finding \og p (J r ) boils down to calculating 

t n dfi(g ■ t) = [ (g- 1 ■ t) n d^{t), 



which is the nth moment of /x at g7L p . This data can be efficiently computed in time polynomial 
in the number of p-adic digits of precision, thanks to the methods of Darmon and Pollack (sec for 
instance [DP061 display 23]). Finally, one recovers the original period via the formula 

J T =p v ^ -C-exppOogpJr), 

where £ is the Teichmullcr lift of the unit part modulo p of the multiplicative integral. Note also 
that v p (J T ) is the sum of the valuations of the ao appearing in the decomposition (|3.7j) . In order 
to find the power series in (|3.8[) , we calculate: 

at+b 

gt-r 2 = g+d _ T2 

(a - cr 2 )t + (b - dr 2 ) 
(a — cT\)t + (b — drx) 

_ b-dr 2 1 + ^f%zf 

= ^p- ( i + jt(-iT(gnrA 

b-dn y gn ^ J 
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where g is the matrix 

'0 -1 
1 / 

Note that for any two matrices g, h we have hg = ~gh~ l . Also, note that for any choice of h, a 
decomposition 

P 1 (Qp) = U^" Z p 



gives rise to another decomposition 



|J(%) ■ Zp- 



Therefore by choosing an appropriate h £ GL2(Q P ) n Mi^K) we can assume that v v {j\ — a) = for 
all a = 0, 1, . . . ,p — 1 and that v p (t2) > 0. 

Note that in order to obtain a power series as in (|3.8|) , the matrices 5 = ( " j ) that we consider 
should satisfy 

(3.9) v p ( > 1. 



d,T2 

The conditions on n and T2 imply that the matrix corresponding to the contribution of P 1 (Q p ) \ Z p 
satisfies p.9[) . and we concentrate on the integral on Z p . Let r > 1 be the largest integer such 
that T2 is congruent to some integer modulo p r . Let t2 € Z be a representative for the class of T2 
(mod p r ). Write also for the representative of ti (mod p l ) in the range 0, . . . ,p l — 1. We can 
then use the decomposition given by the matrices g in the set G = U^^Gj, where 

fi = {(V "" + *") 1 ,-1}, a r+1 = {("; + ' ^)|»_o,...,,-i 

This decomposition consists of exactly p + 1 + r(p — 1) opens. 

3.4. Numerical computations. To test our methods we have written a Sage implementation of 
the above algorithms, modifying an existing implementation written by Robert Pollack which in 
turn adapted part of the code originally written in Magma by Darmon and Pollack ( |DP06j ). The 
code can be found in the second author's webpage. 

Given an elliptic curve of conductor N = pAI and a quadratic field K, the code finds all the 
optimal embeddings of level N of K into M2(Z[i]), and computes the Stark-Heegner period corre- 
sponding to the fixed point of K acting on % p via each embedding to a prescribed precision. The 
Tate parametrization yields the coordinates of the Stark-Heegner point on E(K p ) which are then 
recognized as algebraic coordinates using standard routines. 

Apart from gathering numerical evidence in support of Darmon's conjecture, it is also worth 
remarking that the relative large height of the points thus found would make it impossible to find 
them using naive point search methods. This is, therefore, the only known method to finding points 
of infinite order on such curves. 

The rest of this subsection contains the evidence that we have collected in support of Conjec- 
ture 13.11 Although we do not intend to be exhaustive, we provide examples of curves of small 
composite conductor which satisfy Assumption 13.21 and for which we have been able to positively 
test the conjecture. For each of these curves we consider all the real quadratic fields K of discrim- 
inant D < 200 allowed by the splitting conditions on p and M; for each such field, we consider 
t £ Hp such that H T equals the Hilbert class field of K, and we are able to recognize in all the cases 
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P T as an algebraic point defined over the Hilbert class field of K. For those fields with nontrivial 
class group, we give the relative minimal polynomial feo of the X-coordinate of the point. 



D 


h 


P+ 


13 


1 


(-VT3+l,2\/l3-4) 


28 


1 


(-15\/7 + 43, 150V7 - 402) 


37 


1 


(_|V37+§,i^37-|2) 


73 


1 


V 32 v '° ^ 32' 128 v '° 128/ 


88 


1 


(-¥>!! V22 + D 


97 


1 


f_^5_,/q7 + 123 .375. ./q7_ 4749^ 
V 121 v 31 ' 121 ' 2662 v Jl 2662 / 


133 


1 


/103 92 , /Too 56 N 
\ 9 ' 27 V 100 9 / 


172 


1 


( 1923 11781 /To , 121 \ 
V 1681 ' 68921 v ^° ' 1681/ 


193 


1 


/1885 /Tno , 25885 292175 /TTTo i 4056815^ 
V 288 V ' 288 ' 3456 V 1 3456 / 



Table 1. Points on elliptic curve 15A1, with p = 5 



D 



29 

44 

■53 

92 

137 

149 

197 




209961032 



D 



65 



h D (x) 



J2 , f 61851 
''' ' V 6241 



65 



491926 \ 
6241 / x 



403782 
6241 



65 + 



3256777 
6241 



Table 2. Points on elliptic curve 21 Al, with p = 3 



DARMON POINTS WITH HIGHER CONDUCTOR 



12 




290996167 '' 
7888624 i 



58223 ^ 



21609 ' 3176523 
33533349261 . /-[Qg 



43+ 25921) 
^° ' 21025/ 

46 29744 ^ 



^678412148664608 

147778957920931299317 
12494688311813553741184 




21609/ 
1048806825770477 



678412148664608 ' 

i qo . 30862934493092416035541 
ISO -f 12494688311813553741184 



h D (x) 



16819 



+ - 




83168215321 



83168215321 



66553066916820 '! 
83168215321 / 



Table 3. Points on elliptic curve 33 Al, with p = 11 



DARMON POINTS WITH HIGHER CONDUCTOR 



13 



12565. 



31879 4020800 



SV6- 



. 19321 V u -r lg 32i j 26856 

/IT + 449,2100^41 



12075417 i 



' 7444913385 
, 279945122 

/ 63742245 



2685619 J 

13443) 

)05 . /gj | 14113780406002997 > 



(70 



279945122 



,2805132'^ 



69 + 



526671623 



3312030738382 
1810980922695 



69 + 



3312030738382 
16959961136217 \ 



280513298 ' 3322118988214 v ' 3322118988214 
4398502037370 . /Tg , 19299436937929 



^ 1404725114521 

28304052715333334100 
1664895657706548931 



19 



1404725114521 ' 
121810800584629037164 



'981 3563 /on_ 
k 100 ' 1000 v ° 

7505 310345 



166489565770654893 
1 \ 



64) 



10404' 1061208 v 
oT , 12769 35700 
' 1681 ' 68921 



I 210 
V 1681 

( 9526581863470 

^ 129638878212649 

227155723851142702700 
1476056210913547737643 

/ 4166720 7889580565 



1678197 ^ 
68921 t 




154639065911401 



31843449' 359385165414 



129638878212649 ' 

i on . 5024283358306642389249 
i ~ 1476056210913547737643 

v/l8l-i) 



' 132755895957027 
v 12703756878289 



2 _ 87841 , 85397 
x 9522 ' 6348 

/o7 _ 805260717153160 ^ _ 4207164401474475 . 
' °^ 12703756878289 / X 12703756878289 



/oT , 24540106232139359 
°^ ' 12703756878289 



Table 4. Points on elliptic curve 35 Al, with p = 7 



DARMON POINTS WITH HIGHER CONDUCTOR 
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P+ 
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(!) - 1) 


53 


1 


(3^53+ f;fV 53 + lor) 


77 


1 


/ 5559 /yy I 78911 2040153 /yy , 17804737 \ 
V 55778 V 1 55778 ' 9314926 V 1 9314926 / 


89 


1 


/ 793511 150079871 /on 1 ^ 
V 2401 ' 235298 v oa 27 


101 


1 


( 656788148124048 /TnT , 108663526315570777 
^ 108395925566683225 v iUi 1 108395925566683225' 

432742605985104670344096 /TnT 71551860216079551941383354 \ 
35687772118459783422252125 V 35687772118459783422252125 J 


137 


1 


(83 193 /Toy 1\ 
V81' 1458 v 10 1 27 


149 


1 


/ 41662615293 /T7n 1 802189306199 
I 110013332450 VJ 1 110013332450' 

39791672228037249 /TTn 635290450369692907 \ 
25801976926160750 V J 25801976926160750 J 


152 


1 


( 1915814571 PTo | 24731592007 
I 20670100441 v 1 20670100441 ' 

577303899566856 /qc 7167395643538198 A 
2971761010503011 v °^ 2971761010503011 1 


161 


1 


('62146167667 8395974419456303 /TFT I) 
\ 49710362300 ' 53153799096521000 v u 2) 


188 


1 


( 3178296211866 /Ty 1 22525829850817 
^1135825194001 V ' 1 1135825194001' 






21864116230891316004 /Ty , 148356498531472446055 \ 
1210506836331759751 * ' 1 1210506836331759751 J 






D 


h 


h D {x) 


104 


2 


2 , / 992302702743 ATg 57132410901980 \ 4968445297101 ATfi 1 61480175149213 
^ V 1960400420449 * " U 1960400420449 / 1960400420449 * " U 1 1960400420449 


140 


2 


2 7073157 , 398237221 
x 13924 1 55696 


185 


2 


2 | / 908505900 , /Toe 54207252962 \ 787814100 /Tor , 45005684581 
^ + l 7532677681* 7532677681/ 7532677681 V 1 7532677681 



Table 5. Points on elliptic curve 51 Al, with p = 3 



D 


/i 


P+ 


29 
44 
149 


1 
1 
1 


2-(|V29+f ,|V29+1M) 

V36' 54 v x± 72^ 
/ 41297 ,/TJn , 554429 28371039 /TTn i 340434623 \ 
V 48050 V 1 48050 ' 7447750 7447750 / 



Table 6. Points on elliptic curve 105A1, with p = 3 



DARMON POINTS WITH HIGHER CONDUCTOR 



15 



4. Computation of ATR Darmon points on curves of non-trivial conductor 

In Section |3] we have seen that the algorithm of Theorem 12. 31 can be used in the computation of 
the semi-indchnitc integrals entering the definition of p-adic Darmon points. It is a substitute for 
the continued fractions trick of [DG02] and |DP06j . 

There is another type of Darmon points, called ATR, whose definition also relies in certain semi- 
indefinite integrals. Although the framework is different (e.g., they are points on elliptic curves 
over number fields, and the integrals are complex instead of p-adic), the formal properties satisfied 
by the semi-indefinite integrals are the same in both settings. In the ATR case, the continued 
fraction algorithm over number fields had been used for computing ATR points on curves with 
trivial conductor (cf. |DL03j . [GM12j ). Using a method analogous to that of Section 15721 Theorem 
12.31 also allows for the computation of ATR points on curves with non-trivial conductor. 

To be more precise, let F be a real quadratic number field of narrow class number 1 and let 
O denote its ring of integers. Let E be an elliptic curve over F of conductor 9t, and let L be 
the congruence subgroup consisting of matrices in SL2(0) that are upper triangular modulo 9T. 
Assuming that E is modular, there is a Hilbcrt modular form / of parallel weight two and level 
whose L-series coincides with that of E. Let ujf denote the corresponding T-invariant differential 
2-form onHxW (with L acting on it via the product of the two embeddings of F into R) . 

The following is analogous to Assumption 13.21 

Assumption 4.1. There exists an ideal £> | 9T such that / has eigenvalue 1 for the Atkin-Lehner 
operator W®. 

Let L be the subgroup of PGL2(F) generated by L and the Atkin-Lehner matrix corresponding 
to W$). The previous assumption guarantees that 0, oo S F 1 (F) are f -related. 

Let K be a quadratic ATR extension of F; i.e., a quadratic extension of F that has exactly one 
non-real archimedean place. Suppose that all primes dividing 9t are split in K. We refer the reader 
to |DL03j and |Dar04[ §8] for the definition of the semi-indefinite integrals in this setting. We will 
just mention that they are expressions of the form 

pT py 

J J u+e C/A f , 

where r £ H, x,y £ P 1 (F) are in the same T-orbit, A/ is a certain period lattice that depends on 
/, and ui^ is a non-holomorphic differential easily related to uif. They satisfy analogous properties 
to those of p-adic semi-indefinite integrals; namely 

(i) r ^ = r i: ^ ^ a n 7 e f , r , » G fx, 

(n) r i: + r /; - r i: for au T e «, r , , e fx, 

(in) r s: - r s: - q i: ^ a u n^en, r, S £ f x . 

ATR Darmon points are given by expressions of the form P T = $(A J T J2^°° ^ ), where j T 
belongs to L, $ is the complex uniformization map $ : C/Ae — > E(C), and A is a period conjecturally 
relating and A/. At the cost of replacing P T by a multiple of it, we can assume that j T actually 
belongs to Li(9T) (the notation is as in (|2.3p . with S equal to the archimedean places of F). Then 
the algorithm of Theorem 12.31 computes a decomposition of 7 T of the form 



(4.1) 



7 T = UiLiU 2 L 2 U 3 , 
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where the matrices [/, and Li are of the form 

Ui = ( _ | for some € 0, L; = f ? | for some % € 9T. 

In particular, they belong to T and the same expressions of (|3.4j) and (|3.5j) (changing the multi- 
plicative by the additive notation) express J T ujj in terms of usual double integrals of the 
form ui^ , which in principle can be evaluated by integrating the Fourier expansion of o;^ . 

It is worth remarking that, although the method described above certainly expresses the semi- 
indefinite integrals in terms of definite ones, the running time to directly compute the resulting 
double integrals to a useful accuracy often turns out to be too high. The problem is that if the 
limits of integration are too close to the real axis, then the number of Fourier coefficients needed 
to sum the series to an accurate precision is too high. 

These kind of computational difficulties seem to be inherent to the ATR setting, as they were also 
to some extent present in the initial work of Darmon and Logan |DL03j . In [GM12| some methods 
for accelerating the computation of double integrals in the trivial level case were introduced. The 
authors believe that similar techniques should be applied to the non-trivial level setting in order to 
perform a systematical calculation similar to the one of Section 13.41 

In spite of this, in some simple examples it is possible to directly compute the integrals provided 
by Theorem 12.31 and hence to compute approximations to the ATR points. The following is an 
example of this, which we detail because it provides numerical evidence of the validity of Darmon's 
conjecture in elliptic curves of non-prime conductor. 

Example 4.2. Let F = Q(\/5) and let E be the curve 

y 2 + xy + coy = x 3 - (oj + l)x 2 - (30a; + 45)x - (111a; + 117), oj = — ■ 

The conductor of E is *Xl = (v5 + 6), which has norm 31. This curve was previously considered in 
|Grell| and jGarllj (but note the typo in the displayed equation in both references) . Let a = 1 — y/E 
and let K = F(y/a). The embedding (p: K M^-F) sending oj to the matrix 



W 



3- u> -1 
3-3w -3 + w 



is an optimal embedding of level 91. Under the embedding F R sending y/E to the positive 
square root of 5, the fixed point of W acting on C x is 

r = 0.439291418991 + i ■ 0.353408129753, 

and the image of the unit (—3 + 2uj)oj + 4 — 3oj G 0\ under (p is 



Or 



-4 + 3w 2 - 2lo 
-22 + IQoj 12 - 9a; 



The ATR point attached to (the maximal order of) K is 



J T = I I Ul 



+ 
f ■ 



The determinant of j T is u> + 1, which is a unit. However its upper left entry is not congruent to 1 
modulo 9t. If we let u = — V5 — 2, which is a fundamental unit of F, then —4 + 3a; = u (mod DT). 
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This implies that the matrix 7^. = ("q 1 °) 7t has its upper left entry congruent to 1 modulo 9t. 
We can work with 7^. because 

(4.2) J T = / / w+. 

Observe that det(7^.) = w + 1, which is a unit. The algorithm of Theorem 12.31 works for invcrtiblc 
matrices, not just determinant 1 matrices, and it gives the following decomposition of j' T : 

(ll-w\( 1 W 1 46368+75025u W 1 W 1 -37268-60300^ \ 

V0 1 ) \ 118739-73384^ 1 ) \0 1 )\ -5431444+^3356817 1 )\ 1+w )• 

We use this decomposition to transform (|4.2p into a sum of usual double integrals. The resulting 
integrals have limits not too close to the real axis (the smallest imaginary part is ~ 0.011). Inte- 
grating the Fourier scries with coefficients a m with norm of m up to 180, 000 gives J T to an accuracy 
of approximately 12 digits: 

J T ~ -4.828954817077+ i ■ 4.534696532333. 

There is a point P of infinite order in E(K) having x-coordinate equal to 18883/2420a- 16127/2420 
(this was found using naive search algorithms); let J denote its corresponding image in C/Ae. Then 
the equality 

J T = -13 (mod A B ) 

holds up to the computed accuracy, giving numerical evidence of the equality P T = —2P and, 
therefore, of the rationality of P T . 
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